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ABSTRACT 

A remarkable result of the Wilkinson Microwave Anisotropy Probe (WMAP) observations is that the universe 
was significantly reionized at large redshifts. The standard explanation is that massive stars formed early and 
reionized the universe around redshift z ~ 17. Here we explore an alternative possibility, in which the universe 
was reionized in two steps. An early boost of reionization is provided by a decaying sterile neutrino, whose decay 
products, relativistic electrons, result in partial ionization of the smooth gas. We demonstrate that a neutrino with 
a mass of m„ ~ 200 MeV and a decay time of t ~ 4 x 10 I5 s can account for the electron scattering optical depth 
t w 0. 16 measured by WMAP without violating existing astrophysical limits on the cosmic microwave and gamma 
ray backgrounds. Reionization is then completed by subsequent star formation at lower redshifts. This scenario 
alleviates constraints on structure formation models with reduced small-scale power, such as those with a running 
or tilted scalar index, or warm dark matter models. 

Subject headings: cosmology: theory — intergalactic medium — neutrinos 



1. INTRODUCTION 

The recent observations by the Wilkinson Microwave 
Anisotropy Probe (WMAP) have clarified many issues related 
to the energy content of the universe, and determined key cos- 
mological parameters to unprecedented precision (Spergel et al. 
2003). One of the interesting results of WMAP comes from the 
cross-correlation between the temperature and E-mode polar- 
ization angular power spectra, which has shown that the elec- 
tron scattering optical depth is large, r e = 0.16 ±0.04 (Kogut et 
al. 2003). This high value can be explained by an early phase of 
star-formation, reionizing the universe in a single step at red- 
shift z=17±0.3. A series of papers following the WMAP 
results have described this scenario in detail (Wyithe & Loeb 
2003b; Haiman & Holder 2003; Ciardi et al. 2003; Somerville 
& Livio 2003; Fukugita & Kawasaki 2003; Sokasian et al. 
2003; Cen 2003b; Chiu et al. 2003). The implications of the 
WMAP results for the stellar reionization scenario has been re- 
viewed recently in Haiman (2003). 

In structure formation models in the currently favoured spa- 
tially flat cosmologies with dark energy and cold dark matter 
(ACDM), the earliest large dark matter structures appear at red- 
shift z ~ 20. Stellar reionization is therefore quite natural, pro- 
vided that these early structures are associated with sites of ef- 
ficient star formation and ionizing photon production. Never- 
theless, an optical depth as high as r e = 0.16 requires efficien- 
cies that are higher than expected from the known population 
of galaxies at lower redshifts. Most pre- WMAP studies of cos- 
mological reionization favoured a reionization epoch around 
Z = 7-12 (Haiman & Loeb 1997, 1998; Gnedin & Ostriker 
1997), although significant reionization at higher redshifts was 
predicted in models that included an early generation of mas- 
sive, metal-free stars (Cen 2003a; Wyithe & Loeb 2003a). 

Previous evidence based on the spectra of distant quasars 
have shown that the intergalactic medium (IGM) is highly ion- 
ized at least out to redshifts z ~ 6. However, the most natural 
interpretation of the strong HI absorption seen in the spectra 
of the highest redshift quasars (Becker et al. 2001; Fan et al. 
2003) in the Sloan Digital Sky Survey (SDSS) is that reion- 



ization is being completed near the redshifts of these sources, 
at z ~ 6-7. The thermal history of the IGM contains addi- 
tional evidence that reionization was completed at z < 10 (Hui 
& Haiman 2003; Theuns et al. 2002). 

Post-WMAP analyses (e.g. Haiman & Holder 2003; Cen 
2003b; Wyithe & Loeb 2003b) indicate that such an early gen- 
eration of metal-free stars, or an equivalent 'boost' by a fac- 
tor of > 20 in the efficiency parameters of the z ~ 20 ionizing 
source population, is required relative to that of the z ~ 6 source 
population, in order to simultaneously match the WMAP value 
of T e = 0.16 and the SDSS quasar data. It is worth noting that 
stellar reionization would be much less realistic in cosmolo- 
gies with low small-scale power, in which the onset of struc- 
ture formation is delayed, such as those with a running or tilted 
scalar index (Hannestad et al. 2002; Haiman & Holder 2003; 
Somerville et al. 2003), or warm dark matter models (Barkana, 
Haiman & Ostriker 2001; Yoshida et al. 2003; Somerville et al. 
2003b; Spergel et al. 2003). 

In summary, in order to explain the large value of r e = 0.16, 
the standard stellar reionization scenario requires 1) a hitherto 
unobserved early stellar population, and 2) small-scale dark 
matter structure to be in place by z ~ 20. While these require- 
ments are reasonable, whether they are satisfied remains to be 
tested with future data. In this paper, we pose the question: are 
there alternative possibilities to significantly reionize the uni- 
verse at z^ 20? One such alternative that avoids the above re- 
quirements is reionization by decaying particles. Sciama (1982) 
suggested the possibility of photo-reionization from radiatively 
decaying particles where the decaying particle has a mass of 
few tens of eV. Related issues were discussed thoroughly in 
numerous papers (early studies include Cabibbo et al. 1981; 
Salati and Wallet 1981; Rephaeli and Szalay 1981). However, 
these types of particle physics candidates, such as light photi- 
nos and gravitinos, are now excluded by modern accelerator 
data (Hooper and Plehn 2002; Dreiner et al. 2003). Another 
interesting candidate was a decaying active neutrino (Sciama 
1990) which was proposed both as a dark matter candidate 
and also to reionize the universe. However, current cosmolog- 
ical data excludes such massive neutrinos (Adams et al. 1998; 
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Bowyer et al. 1999; Spergel et al. 2003). 

Let us imagine the general case of a cosmologically signifi- 
cant particle, with the property that it decays around a redshift 
z <~ 20. The decay products could have energies high enough 
to reionize the light elements. If this hypothetical particle is 
sufficiently abundant, then it may explain the high value of the 
optical depth, without the need for stars to reionize the universe 
early. In such a scenario, normal stellar populations could form 
at lower redshifts (z < 10), with the usual efficiencies, seeded in 
dark matter structures on the better understood scales of dwarf 
galaxies, and account for the completion of reionization epoch 
that appears to be occurring at z < 10. 

In this paper, we turn to particle physics and describe a spe- 
cific example for such a particle that is both viable (on par- 
ticle physics grounds) and satisfies existing astronomical con- 
straints. In practice, the choice of a specific particle candidate 
turns out to be severely limited by astronomical observations, 
such as big bang nucleosynthesis (BBN), the diffuse gamma- 
ray background (DGB), and the stringent limit on the deviation 
of the spectrum of the cosmic microwave background (CMB) 
from a pure Planck shape. Here we present a model which re- 
lies on a minor extension of the standard model: introducing a 
heavy sterile neutrino which mixes with one of the active neu- 
trinos. Heavy sterile neutrinos are expected in extensions of the 
standard model, in order to give masses to the active (light) neu- 
trinos through the see-saw mechanism (Yanagida 1980; Mo- 
hapatra and Senjanovic 1980). The neutrino masses we will 
consider below are somewhat smaller than invoked for the see- 
saw mechanism. Our particle should therefore be regarded as a 
new independent sterile neutrino (see e.g. Hansen 2001 ; Dolgov 
2003, for reviews on heavy sterile neutrinos in cosmology). For 
the purpose of reionization, the main novelty of this particle is 
that its decay products are energetic electrons, rather than pho- 
tons. As we shall see below, this makes it possible to achieve 
reionization without violating other astrophysical observations. 

The rest of this paper is organized as follows. In § 2, we de- 
scribe the candidate particle we propose, quantify the expected 
amount of reionization, and show that the particle obeys various 
astrophysical constraints. In § 3, we discuss the cosmological 
implications of this new particle, and in § 4, we offer our con- 
clusions. Throughout this paper, we adopt the background cos- 
mological parameters as measured by the WMAP experiment, 
n m = 0.29, il A = 0.71, Q b = 0.047, h = 0.72 (Spergel et al. 2003, 
Tables 1 and 2). 

2. DECAYING STERILE NEUTRINO 
2.1. Production and Decay Mechanisms 

We will now discuss some details of the production and de- 
cay mechanism of the heavy sterile neutrino we propose. Our 
main goal is to demonstrate that particle physics can provide a 
viable decaying particle whose decay products (energetic elec- 
trons in our case) can reionize the universe early, without vio- 
lating other observations. We do not attempt a comprehensive 
analysis of all possible particle physics candidates here; indeed, 
we expect that there could be other particle candidates which 
could similarly account for early reionization. 

Let us consider a standard two-neutrino oscillation scheme, 
where one of the active neutrinos, v a = v e , or v T , mixes with 
a heavy mainly sterile neutrino, v s , 



v a = cos9 v\ + sin# V2 , 
v s = - sin 9 v\ + cos 9 i>2 , 



where v\ and are assumed to be the light and heavy mass 
eigenstates respectively, and 9 is the vacuum mixing angle. Let 
the light neutrino be massless and the heavy neutrino have mass 
m. We will consider very small mixing angles, and hence refer 
to the light neutrino mass eigenstate as the active neutrino and 
the heavy one as the sterile neutrino. 

If the sterile neutrino has a mass in the range 140-500 MeV, 
then the dominant decay channel (Astier et al. 2001) is 



l + TT, 



(2) 



where the lepton, can be either an electron, a positron, or a 
neutrino. For larger masses, new decay channels open up (in- 
cluding kaons, eta, etc.), and the considerations below become 
more complicated. For smaller masses the decay products in- 
clude too few photons and electrons for our purpose. The rate 
for the reaction 7ro — ► vv was calculated by Fischbach et al. 
(1977) (see also Kalogeropoulos et al. 1979), and the result was 
converted to the reaction v s — > v a + ttq by Dolgov et al. (2000). 
Here we also need to include the charge current reactions 



e +7T+, 
e + + 7r_, 



(3) 



and we find the lifetime 
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(4) 



where r 2 w 0.2, Gf is the Fermi constant, = 131 MeV, and 
the masses M = m/m 7T are all in units of the charged pion mass, 
mn = 139.7 MeV. When we express the decay time, r, in units 
10 15 sec « 30 Myr, then this equation becomes 



t 15 M(M 2 - l)sin 2 <9 = 1.2 x 10 



-24 



(5) 



The reactions (3) directly provide us with a significant number 
of highly relativistic electrons and positrons, which we can use 
to reionize the universe early. Some of the energetic positrons 
could annihilate with the thermal (and also with the much more 
dilute, energetic) electrons, but the resulting photons would 
have mean free path larger than the Hubble distance. Such anni- 
hilation would change our result slightly at worst, by reducing 
the number of ionizations. 

Let us now turn to the production and abundance of the ster- 
ile neutrinos. These neutrinos are produced through standard 
neutrino oscillations while they are still relativistic, at a tem- 
perature of about 7GeV, when there are about 90 degrees of 
freedom in the plasma (Kolb and Turner 1990). For mixing 
with electron neutrinos, the fraction of sterile neutrinos to ac- 
tive ones is (Dolgov & Hansen 2001) 



— =4.0x lO 9 sin 2 0M. 



(6) 



(1) 



For mixing with muon or tau neutrinos, the numerical coeffi- 
cient in front is instead 5.6. If the sterile state is mixed with 
both electron and muonic (tau) neutrinos, then the produced 
number increases correspondingly. We thus assume that there 
is no thermal relic abundance of these sterile neutrinos. Spe- 
cific models with such properties exist (Bento & Berezhiana 
2001; Berezhiani & Mohapatra 1995), however, other produc- 
tion mechanisms can produce a significantly higher abundance 
of sterile neutrinos (Sarkar 1996; Mohapatra and Sciama 1998). 
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2.2. Reionization Efficiency 

Let us next estimate the amount of hydrogen which can be 
reionized by the decay products. The universe is transparent to 
the vast majority of photons from the process no — ► 27, since 
these photons have energies above m v /2. Since the abundance 
of these sterile neutrinos is about 15 orders of magnitude below 
that of a thermalized neutrino, one easily sees (Kolb and Turner 
1990) that the diffuse gamma background is many orders of 
magnitude larger than the photons resulting from both ttq and 
tt^~ . We can therefore ignore this fraction of the decay products 
and consider the fate of the relativistic decay electrons instead. 
These electrons are injected into the IGM at redshift z ~ 20 (i.e. 
at the redshift where the decay time equals the age of the uni- 
verse) with energies of about E e = (m^jT)(M- 1), i.e. between 
< E e < 180 MeV. The relevant interactions for such electrons 
with the background gas and with the CMB have been summa- 
rized in, e.g. Oh (2001). For our purposes, the two important 
processes are inverse Compton scattering with CMB photons, 
and collisional ionization of neutral hydrogen. The mean free 
path to Compton scattering is 

\c = (<y T n i r x ~4x 10 17 [(l+z)/21]- 3 cm, (7) 

where a T is the Thomson cross section and n 1 is the CMB pho- 
ton number density. The mean free path to collisionally ioniz- 
ing hydrogen is 

Ah = {o H n H T X ~ 1-4 x 10 24 (£ c /100MeV)[(l + z)/21]- 3 cm, 

(8) 

where 07/ is the electron impact ionization cross section (Shull 
& van Steenberg 1985, equation A2; we have ignored a small 
logarithmic correction), and nj = 2x 10~ 3 [(1 + z)/21] 3 cm~ 3 is 
the neutral hydrogen density. Since Ac -C Xh, the electrons will 
undergo repeated Compton scatterings with the CMB photons. 
At each scattering event, the CMB photon will be up-scattered 
to the energy E 1 ~ (4/3) 7 2 £ 7 , ~ 300(£ e /100MeV) 2 [(l + 
z)/21] eV, where < 7 < 360 is the Lorentz factor for the 
electron, and £ 7 ,o = 5 x 10~ 3 [(1 +z)/21] eV is the mean (undis- 
turbed) CMB photon energy. The up-scattering of the CMB 
photons is the first step (of 3 steps) towards reionization. The 
electron will loose the same amount of energy E 7 , and will 
cool after N ~ E e /E 1 = 3 x 10 5 (£ e /100MeV)- 1 [(l +z)/2l]~ l 
scatterings, corresponding to the cooling time 17 ~ NXc /c ~ 
6 x 10 4 (£ , e /100MeV)- 1 [(l +z)/21]" 4 yr. The cooling time for 
the electron due to collisionally ionizing hydrogen, on the other 
hand, is t h - 2.9 x 10 9 (£ e /100MeV)- 1 [(l +z)/21]" 3 yr. 

By equating tq = th we infer that decay electrons whose en- 
ergies are above E e > 1[(1 +z)/21] -1 ' 2 MeV will loose most 
of their energy to up-scattering the CMB photons. We con- 
clude that only electrons that are produced by sterile neutri- 
nos very near (within ~ 2 MeV) the pion mass could spend a 
significant fraction of their initial energy on direct collisional 
H-ionizations. However, CMB photons that are up-scattered 
by electrons with energies 20[(1 +z)/21] _1 / 2 MeV < E e < 180 
MeV will have final energies 

13.6eV<£ 7 <900[(l+z)/21]eV. (9) 

These photons then spend a significant fraction of their en- 
ergy on reionizing hydrogen. At the low (UV) end of this 
energy range, the Comptonized photons directly photoionize 
the neutral hydrogen atoms. Despite the steep decrease of 
the photo-ionization cross-section (a v oc i/" 3 ), even at the high 
(soft X-ray) end of the range, the mean free path of the pho- 
tons to photo-ionization is a small fraction (cov«#fHubbie) _1 ~ 



0.2(£ c /180MeV) 3 [(l +z)/21] 3 / 2 of the Hubble distance. As a 
result, each photon ionizes a hydrogen atom, and creates a fast 
secondary photoelectron. The photo-ionization and production 
of the energetic photoelectron is the second step towards reion- 
ization. At these low (non-relativistic) electron energies, in- 
verse Compton scattering is negligible compared to collisional 
ionization, and a fraction <~ 1 /3 of the initial energy of the pho- 
toelectron is subsequently spent on ionizing hydrogen (Shull & 
van Steenberg 1985; Oh 2001; Venkatesan et al. 2001). This 
collisional ionization is the third step towards reionization. 

We conclude that a significant fraction (> 1/3) of the de- 
cay electrons can produce hydrogen ionizations. The maximal 
number of hydrogen atoms that could be ionized by a single 
electron, on energetic grounds alone, would be \ = E e /l3.6eV. 
This can be written as % = 5.1 x 10 6 (M- 1), where M is again in 
units of the pion mass. Thus the number of reionized hydrogen 
atoms is given by 

(10) 



prod 



prod 



where we have introduced an efficiency factor, e < 1 , and n e 
is the number density of electrons (plus positrons) from reac- 
tion (3). Using «i/n 7 = 6 x 10" 10 and n a = 3/ll« 7 , we find 
that the required efficiency is 

e = 0.10r 15 (M+l) . (11) 
For mixing with muon (or tau) neutrinos the factor in front 
is 0.071. Note that this efficiency is quite reasonable: using 
e = 1/3 (see above discussion) and assuming mixing with an 
electronic neutrino, one can reionize half of the baryons with 
r 15 (M+ 1) = 6.7 (see Fig. 1). 
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FIG. 1 . — Parameter space (decay time vs. mass) for a decaying sterile neu- 
trino mixed with a electron neutrino, using a reionization efficiency e= 1/3. 
The right axis shows the redshift corresponding to the given decay time. The 
solid curve corresponds to ionizing half of the baryons, while the dashed and 
dotted curves correspond to ionizing a smaller fraction of 1/3 and 1/5, respec- 
tively. The star indicates the life-time corresponding to the redshift where the 
universe should be reionized to the fraction 33% in order to give an optical 
depth of 0.16. 

2.3. Astrophysical Constraints 

Let us next turn toward the observational consequences of the 
decaying sterile neutrinos. First, we can consider how many 
sterile neutrinos were produced, and whether they will affect 
the production of light elements (BBN), photon decoupling, or 
supernova cooling. Using equations (5) and (6), we find 

IK ;5 xl0- 15 . (12) 



r 15 (M2-l) • 
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For masses in the range considered (140 MeV< m < 500 MeV) 
and for a lifetime several orders around r = 10 15 sec, this is 
such a small number that none of the three processes above will 
be affected in any way (Kainulainen et al. 1991; Sarkar 1996; 
Dolgov et al. 2000; Dolgov 2003). The neutral pions decay im- 
mediately into two energetic photons; however, the observed 
diffuse gamma ray background is several orders of magnitude 
higher than the signal from the decay products (Kolb and Turner 
1990). Furthermore, the mixing angle is so small (see eq. 5), 
that no terrestrial experiment can detect the oscillations. 

A potential observable from reionization at high redshift is 
the spectral distortion of the CMB. The CMB photons have 
been observed by the COBE satellite to obey a nearly perfect 
Planck distribution. The COBE results translate into an up- 
per bound on the energy deposited by the electrons, usually 
expressed in terms of the limits | y |< 1.5 x 10" 5 and | p |< 
9 x 10~ 5 on the Compton y parameter and the dimensionless 
chemical potential, respectively (Fixsen & Mather 2002; Fixsen 
et al. 1996). As we showed above, in our scenario, a significant 
fraction of the total energy of the decay electrons is transfered 
to the CMB. For a sterile neutrino near the pion mass thresh- 
old (1 < M < 1.02), for which the decay electrons collisionally 
ionize (and heat) the IGM directly, this transfer occurs by the 
non-relativistic Compton coupling between the reionized IGM 
and the CMB. For high energy decay electrons, the transfer oc- 
curs in two steps. CMB photons are first up-scattered to UV 
or soft X-ray energies, effectively removing a small fraction 
of the photons from the COBE range, corresponding to a p- 
distortion. Second, this energy is transfered back to the COBE 
range by non-relativistic coupling between the reheated IGM 
and the CMB (corresponding to a y-distortion). 

We therefore compute the total energy in the primary decay 
electrons (plus positrons), p e , and compare this to the total en- 
ergy of the CMB, p 7 . Using the CMB temperature T « 5 x 
10" 3 eV (at redshift z=20) we find p e / Pl « 7 x 10- 6 /(ti 5 (M + 
1)). Depending on the actual fraction (0.1-0.9) of the mean 
electron energy that goes into heating the IGM (rather than go- 
ing into the ionizations and excitations), the number can be 
even smaller. We conclude that both the p- and y-distortions 
are sufficiently small (both < p e j p 7 ) that there is no disagree- 
ment with the current upper limit from COBE. However, the 
decaying neutrino model presented here should predict a char- 
acteristic shape for the distorted CMB spectrum that could be 
detected in future measurements improving the sensitivity by 
nearly two orders of magnitude (Fixsen & Mather 2002). Fi- 
nally, we note that in the stellar reionization case, where the 
ionizing sources are strongly clustered, the angular fluctuations 
of the spectral distortion can place interesting limits on models 
(Oh et al. 2003). Since the electrons from the decaying neutri- 
nos are smoothly distributed in space, our scenario avoids this 
constraint by construction. 

3. DISCUSSION 

Let us briefly summarize the scenario envisioned in this pa- 
per. Sterile neutrinos with masses of about a few hundred MeV 
were created in the early Universe with an abundance about 
five orders of magnitude below the number density of baryons. 
A large fraction, 20-30% of the neutrino's mass results in 
decay-produced energetic electrons and positrons. Depending 
on their energies, these electrons then produce a UV/soft X-ray 
background by inverse Compton scattering CMB photons. The 
newly created energetic photons then ionize the IGM either by 



direct photo-ionization, or by the collisions between the fast 
secondary photoelectrons and the neutral hydrogen atoms. For 
a sterile neutrino mass near (within <~ 2 MeV) the pion mass 
threshold, the ionization is by direct electron-hydrogen colli- 
sions. As a result, any sterile neutrino with a mass between 
140 < M < 500 MeV would produce relativistic electrons that 
spend > 1 /3 of their energy on partially reionizing the universe. 

In this scenario, the sterile neutrinos produce electrons at a 
rate oc exp(-f/r). For sterile neutrino masses near the lower 
mass ranges considered here, which result in electrons with 
E < 100 MeV and ultimately produce UV photons, the ion- 
ization time is always shorter than the Hubble time, includ- 
ing at redshift z ~ 1000. As a result, reionization would 
track the neutrino decay rate, and neutrinos with these masses 
would likely lead to a "tail" of the recombination epoch, rather 
than an episode of 2nd reionization at lower redshift z ~ 30 
(as in Adams et al. 1998). Such a scenario may be in con- 
flict with existing weak constraints from CMB polarization 
anisotropies (Kogut et al. 2003; Bean et al. 2003). In com- 
parison, for more massive neutrinos, reionization would be a 
three-step process, as explained above: (1) the relativistic elec- 
trons produce X-ray photons from CMB; (2) the X-ray photons 
photo-ionize hydrogen atoms, and (3) the fast ~keV photoelec- 
tron collisionally ionizes many more hydrogen atoms. Of these, 
(1) and (3) last only a small fraction of the Hubble time for any 
choice of our model parameters. However, for the upper energy 
range we considered for the initial relativistic electron, at red- 
shift z ~ 1000, (2) will last for ~ 10 8 yr, a time comparable to 
the age of the universe at z ~ 30. As a result, for these energetic 
electrons, reionization will occur near z ~ 30, producing a 2nd 
reionization at these lower redshifts. 

We note further that the universe cannot be fully reionized 
at high redshift in this way (Oh 2001; Venkatesan et al. 2001). 
While the IGM is nearly neutral, the high energy electrons will 
spend about 1/3 of their initial energy on hydrogen ionizations, 
and very little energy will be deposited into heating the ther- 
mal electrons (Shull & van Steenberg 1985). Once the ioniza- 
tion fraction rises to 50%, the fraction of the electron energy 
spent on ionizations drops to only few percent, and the fraction 
spent on heating of the thermal electrons rises to 90% (Shull 
& van Steenberg 1985). It is therefore a reasonable guess that 
the electrons will drive the ionization fraction up to x e <~ 0.5. 
Note that at redshifts of 15-20, the recombination time is a 
fraction 0.5-0.3 of the Hubble time; we therefore expect re- 
combinations to allow this level of partial ionization to remain. 
A universe reionized to 30(50)% at redshifts of 7 < z < 31(25), 
but fully reionized at < z < 7, has an optical depth of r e = 0. 16 
(in our adopted ACDM cosmology, conservatively ignoring any 
contribution to the opacity from He+ or He++ which recombine 
more rapidly than hydrogen). We conclude that a scenario in 
which a decaying particle produces an early boost of reioniza- 
tion, followed by full reionization by stars, is therefore plausi- 
ble, and can account for the WMAP optical depth. 

Our scenario implies that warm dark matter (WDM) models 
are not ruled out. This is not in contradiction with the results of 
Barkana et al. (2001), Spergel et al. (2003) and Yoshida et al. 
(2003), since all of these studies assumed that reionization was 
done by stars alone. The possibility of a sterile neutrino saving 
the option of WDM is rather peculiar, because the only (to our 
knowledge) natural WDM candidate is, in fact, another sterile 
neutrino (Dodelson & Widrow 1994; Dolgov & Hansen 2001; 
Abazajian et al. 2001), which is constrained to have a mass of 
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a few keV (Abazajian et al. 2001; Hansen et al. 2002). One 
should mention that more inventive WDM candidates certainly 
exist (see e.g. Baltz & Murayama 2001; Boehm et al. 2002). 

Thus, with a minimal extension of the standard model with 
two sterile neutrinos (one with a mass of about a keV, and an- 
other with a mass of a few hundred MeV), one obtains the 
right abundance of dark matter, and achieves an early boost 
of reionization. WDM models may be attractive in helping 
to solve some possible problems of CDM on small-scales, al- 
though many issues remain to be clarified (see summary in, e.g., 
Haiman et al. 2001). 

A very important question is whether we will be able to dis- 
tinguish between early stellar reionization models, and the sce- 
nario described above. It was shown by Kaplinghat et al. (2003) 
that two years of polarization data from WMAP can distinguish 
between a scenario with full reionization up to z = 6.3, and a 
model with an additional 50% reionization between z = 20 and 
6.3. Planck will do even better in just one year, namely distin- 
guishing between 0% or 6% reionization between z = 20 and 
6.3. For related studies, see e.g. (Doroshkevich et al. 2003; 
Naselsky & Chiang 2003). Specifically, Haiman & Holder 
(2003) and Holder et al. (2003) showed explicit examples of 
different reionization histories which produce the same optical 
depth r = 0.16, but can be distinguished from each other by 
Planck polarization measurements at > 3a significance. Based 
on these studies, we are confident that the model discussed 
above can be constrained or confirmed within a few years. The 
reionization history in the decaying neutrino model could pos- 
sibly be mimicked by ionization caused by an early generation 
of stars that produces X-rays and therefore fast photoelectrons 
(Oh 2001; Venkatesan et al. 2001). However, this "degener- 
acy" can be broken by examining small angular scale Sunyaev- 
Zeldovich fluctuations, which should be strong in the stellar 
case (Oh et al. 2003), but nonexistent in our scenario. 



There are many issues we have not touched upon in this sim- 
plified treatment. For example, hydrogen recombinations are 
significant at high redshifts, and one would have to convolve 
the recombination rate in the IGM with the neutrino decay rate, 
in order to compute the actual ionization fraction as a func- 
tion of redshift. While we do not expect such a computation 
to change our main conclusions, it will allow predictions to be 
made about the precise evolution of the ionization and thermal 
history of the IGM, as well as the resulting spectral signatures 
on the CMB. 

4. CONCLUSIONS 

In this paper, we studied the possibility of a reionization his- 
tory that consists of two completely independent steps. The first 
step is provided by fast electrons from a decaying sterile neu- 
trino. The abundance and mass of this sterile neutrino is such 
that up to 50% of the universe can be reionized around a redshift 
of z = 25. Subsequently, standard (population II) star formation 
can complete reionization at z ~ 6-7. This scenario can give 
an electron scattering optical depth of r = 0.16, in agreement 
with WMAP observations. In such a scenario, the possibility of 
warm dark matter still exists, since stars only need to complete 
the reionization around redshift of z = 7 to comply with the ob- 
served Gunn-Peterson opacities of high redshift (z ~ 6) quasars. 
The model proposed here can be tested in future studies of the 
CMB, both by inferring the reionization history (evolution of 
the ionized fraction) from post-WMAP polarization data, and 
by studying the deviation of the CMB spectrum from a pure 
black body shape. 
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